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We present a search for differences in the oscillations of antineutrinos and neutrinos in the Super- 
Kamiokande -I, -II, and -III atmospheric neutrino sample. Under a two-flavor disappearance model 
with separate mixing parameters between neutrinos and antineutrinos, we find no evidence for a 
difference in oscillation parameters. Best fit antineutrino mixing is found to be at (Am 2 , sin 2 2#) = 
(2.0 x 10 _3 eV 2 , 1.0) and is consistent with the overall Super-K measurement. 
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PACS numbers: 14.60.Pq, 96.50.S- 



As the parameters outlining the standard neutrino 
oscillation framework become increasingly well known, 
searches for sub-leading and possibly symmetry-breaking 
effects become possible. If the value of $13 is non-zero, for 
instance, it becomes possible to search for CP-violation 
effects in the neutrino system via differences in the os- 
cillation probabilities of neutrinos and antineutrinos. In 
this paper we consider the possibility that the survival 
probability P(V M — > v^) is governed by a different mass 
splitting or mixing angle compared to P(t' Al — > v^)- This 
is not considered in most oscillation studies as the mass 
splitting and mixing angle are expected to be identical 
for neutrinos and antineutrinos by CPT symmetry. An 
inequality of these probabilities, in the absence of matter 
effects, could signal new physics. For atmospheric muon 
neutrino disappearance, which is predominantly oscilla- 
tion into tau neutrinos, the matter effect is expected to 
be small and the mixing parameters for muon neutrino 
and antineutrino disappearance appearance should be 
the same. The MINOS experiment, which is sensitive to 
neutrino oscillations at the atmospheric scale, and which 
can determine the sign of muons by magnetic bending, 
has observed antineutrino disappearance P at a best- 
fit value of Am 2 nearly 50% larger than previous mea- 
surements made using neutrinos [2H5]. Though not in 
the realm of atmospheric mixing, the MiniBooNE exper- 
iment has similarly observed a discrepancy between its 
neutrino [5] and antineutrino ?i| data. Therefore, further 
tests of differences between neutrinos and antineutrinos 
using atmospheric data are well motivated. 

Super-Kamiokande (Super-K, SK), described below, 
cannot distinguish v from D on an event by event ba- 
sis so potential differences in their oscillations would ap- 
pear in the atmospheric neutrino sample in a statistical 
way. Notably, the neutrino and antineutrino cross sec- 
tions differ by a factor of two to three depending on the 
neutrino energy. The ratio of the v and v atmospheric 
fluxes is similarly energy dependent [5]. For these rea- 
sons, even in the absence of CPT-violating oscillations 
the relative numbers of each species are expected to dif- 
fer. Kinematic considerations can also induce differences 
in the products of neutrinos and antineutrino reactions, 
enhancing the sample purity of one or the other. For 
instance, the absorption of 7r~ on 16 nuclei in water 
tends to enrich the neutrino component of samples that 
are sub-divided based on their number of decay electrons. 
If CPT-violating oscillations are present in the data, the 
distortion of the zenith angle distribution characteristic 
of disappearance would appear at different energies 
and path lengths (different oscillation frequencies) be- 
tween neutrinos and antineutrinos. Since Super-K can 
only observe the total distribution, a potential signal 
would appear as a distortion consistent with the com- 



position of separately oscillated spectra. In this context, 
antineutrinos, which have the smaller cross section, are 
expected to provide a weaker oscillation constraint than 
neutrinos. 

In this Letter we consider ad-hoc CPT-violating oscil- 
lations by testing separate two-neutrino disappearance 
models for neutrinos and antineutrinos: 

2 a fAm 2 L\ 
P(f„ -> = l-sm 2 2#sm( J 

Pfo -> Vy) = 1 - sin 2 2#sin , (1) 

where L is the neutrino path-length and E is the neu- 
trino energy. In the presence of matter, additional 
neutrino-electron scattering induces a CPT-violating-like 
difference between the neutrino and antineutrino survival 
probabilities [§1 [TU], particularly when #13 is non-zero. 
Recent data suggest [TlTfl3] that #13 is small and there- 
fore matter-induced corrections to Eqn. [I] are expected 
to be sub-dominant. Changes to the fit results induced 
by a three-flavor treatment are briefly considered below. 
Although the presence of both matter and solar mixing 
terms is expected to drive — \ v e transitions below 
about 1 GeV even if 6*13 = 0, the oscillation frequency in 
this domain is high enough that the effects are averaged 
out by the detector resolution. Since the dominant at- 
mospheric disappearance effect is seen at higher energies, 
a two-neutrino scheme is the focus of this study. 

Super-Kamiokande is a water Cherenkov detector lo- 
cated in Japan's Gifu prefecture and situated at a depth 
of 2700 meters water equivalent. It is comprised of two 
concentric, optically separated cylinders: an inner de- 
tector (ID) viewing a 22.5 kton fiducial volume and an 
outer detector (OD) used primarily as a veto. During the 
first run of the detector, SK-I, the walls of the ID were 
lined with 11,146 inward- facing 20-inch photomultiplier 
tubes (PMTs). The two subsequent run periods, SK-II 
and SK-III, had 5,182 and 11,129 ID PMTs, respectively 
and the PMTs have been encased in fiber-reinforced plas- 
tic shells. The OD has been instrumented with 1,885 
outward- facing 8-inch PMTs throughout. More details 
on the detector and its calibration may be found in |14) . 

The atmospheric neutrino data are divided into three 
categories. Fully contained (FC) events deposit all of 
their light in the ID, partially contained (PC) events ad- 
ditionally have an exiting particle that deposits energy in 
the OD, and upward-going muon (Up-/i) events are pro- 
duced by neutrino interactions in the rock beneath the 
detector. Up-/x events are required to have a minimum 
pathlength of seven meters and are classified as stopping 
or through-going. Through-going Up-/i events are fur- 
ther sub-divided into "showering" and "non-showering" 
based on [T5] . 
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FIG. 1: SK-I+II+III lepton momentum (first column) and zenith angle distributions of the event samples used in the analysis. 
Black dots represent the data with statistical errors, the histogram is the oscillated MC expectation at the best fit point with 
the shaded region showing the antineutrino composition. 



The present analysis uses data accumulated during the 
first three SK run periods. SK-I (1996 to 2001) FC and 
PC events correspond to 1489 live-days with 1646 days of 
Up-/x livetime. SK-II (2002 to 2005) had 799 FC/PC and 
828 live-days of Up-/i events. The FC and PC livetimes 
during SK-III (2005-2007) were 518 days and that for 
Up-// was 636 days. 

During the analysis, FC and PC events are further 
divided. Fully contained events are separated into sub- 
GeV ( E OTS < 1.33 GeV ) and multi-GeV ( E vis > 1.33 
GeV ). These samples are then separated based on their 
number of reconstructed Cherenkov rings into single- and 
multi-ring topologies. Pattern-identification of single- 
ring events is used to separate them into e-like and /i-like 



categories. This technique is applied to the multi-ring 
sample using the most energetic Cherenkov ring. The 
sub-GeV single-ring e-like and /i-like samples are also di- 
vided based upon their number of decay-electrons. A 
7r°-like sample is also extracted from the single-ring e- 
like events [IS]. Partially contained events are separated 
into "OD stopping" and "OD through-going" categories 
based on the amount of Cherenkov light observed in the 
OD at the exit point. 

Since the physical configuration of the detector and its 
reconstruction performance varies among the SK run pe- 
riods, separate 500 year MC samples are used for each. 
A "pulled" \ 2 \16^ based on a Poisson probability distri- 
bution is used to compare the data against the MC: 
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In this equation n indexes the data bins, E% K% is the MC 
expectation for SK-z, and Of l K% is the number of observed 
events in the n th bin during SK-i. The effect of the i th 
systematic error is introduced via the error parameter 
€i and where the latter is the fractional change in 
the MC expectation of bin n introduced by a 1-sigma 
shift in its systematic error, ui. The data and MC are 
divided into 420 bins for each of the SK run periods when 
computing these systematic errors, but are later merged 
as above to ensure the stability of the fit function against 
sparsely populated bins. In total, 420 bins are used to 
compute the value of \ 2 - 

Equation (12]) is minimized with respect to the ac- 
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cording to = 0, yielding a set of linear equations in ej 
that are solved iteratively [16]. Following this procedure 
a x 2 value is computed for each point in the oscillation 
parameter space. The global minimum x 2 is defined as 
the analysis' best fit point. 



Parameter 


Best Fit 


90% C.L. 


Three-Flavor 


Am 2 (eV 2 ) 
Am 2 (eV 2 ) 
sin 2 26> 
sin 2 26 


2.1 x 10" 3 
2.0 x 1CT 3 

1.0 

1.0 


1.7 , 3.0 ]xl0~ 3 
1.3 , 4.0 ]xl0~ 3 
[ 0.93, 1.0 ] 
[ 0.83, 1.0 ] 


1.7 , 3.3 ]xl0~ 3 
1.2 , 4.0 ]xl0~ 3 
[ 0.93, 1.0 ] 
[ 0.78, 1.0 ] 



TABLE I: Best fit information for the four parameter fit 
to the SK-I+II+III data. The 90% C.L. column represents 
bounds taken from single-parameter Ax 2 distributions in 
which the remaining three parameters have been minimized 
over. The third column shows the 90% C.L. allowed region 
when the effects of #13 and <5 cp are considered in the fit (see 
text). 



The 120 sources of systematic uncertainty considered 
in this analysis are separated into two classes: those that 
are common throughout the SK run periods and those 
that are dependent upon a particular detector geometry. 
Common systematic errors stem from uncertainties in the 
neutrino interaction cross sections, nuclear effects, and 
the atmospheric neutrino flux. Independent systematic 
errors are related to detector performance and include 
uncertainties in the event reconstruction and reduction. 
A complete list of the systematic errors used here is pre- 
sented in [12] . 

Antineutrino oscillations are considered independently 
of neutrino oscillations over a four-dimensional oscilla- 
tion space with two parameters for each: (Am 2 , sin 2 26) 
and (Am 2 , sin 2 26 ) ). All parameters are varied simul- 
taneously on a grid of 50 x 35 points in the antineu- 
trino plane and 20 x 10 points in the neutrino plane. 
The neutrino (antineutrino) parameter space is taken over 
1.0(0.7) x 10~ 3 < Am 2 < 5.0(8.0) x 10~ 3 eV 2 and 
0.85(0.65) < sin 20 < 1.0(1.0), comprising an area en- 
compassing the current allowed values of these param- 
eters [SJ [12]. Minimizing the \ 2 function in Eqn. [2] 
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FIG. 2: Allowed regions for the antineutrino mixing pa- 
rameters for the SK-I+II+III data set. The 68%, 90%, and 
99% allowed region appear in thin, medium, and thick lines 
respectively. The shaded region shows the 90% C.L. al- 
lowed region for antinuetrino disappearance in an antineu- 
trino beam from MINOS pQ. A solid point denotes the 
location of the best fit from that analysis. 



over this parameter space, the best fit is found with 
X 2 = 468.4 for 416 degrees of freedom. Table [i] sum- 
marizes the fit information. Since the oscillations of at- 
mospheric neutrinos are sensitive to the effects of #13 and 
5 cp , an additional analysis has been performed assuming 
distinct three-flavor mixing between neutrinos and an- 
tincutrinos. Though this model has no strong theoretical 
motivation, the potential impact of a three-flavor frame- 
work on the SK allowed regions is provided in the table 
for reference. In both fits no difference between antineu- 
trino and neutrino mixing is found in the data. 

The combined data overlaid with the MC expectation 
at the best fit point are shown in Fig. [T] Antineutrinos 
have been oscillated independently of neutrinos in the 
expectation and are shown as the shaded portion of the 
histogram. Figure [2] shows the allowed regions at several 
C.L. in the antineutrino plane. The neutrino parameters 
have been minimized over and the contours have been 
drawn for a standard A% 2 distribution with two degrees 
of freedom. The 90% C.L. regions from the antineutrino 
and neutrino parameters overlaid with the allowed region 
from the standard two-flavor analysis, where the mixing 
parameters are required to be identical for neutrinos and 
antineutrinos, appear in Fig. [3] The best fit point of 
the standard analysis, (Am 2 , sin 2 26>) = ( 2.1 x 10~ 3 eV 2 , 
1.0 ), is consistent with the best fit from the antineutrino 
separated analysis. 

At each point in the plane of Fig. [2] the best fit to 
the data may lie at a point in the neutrino parameters 
that does not correspond to equal neutrino and antineu- 
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FIG. 3: Allowed regions for atmospheric mixing parameters 
from the SK-I+II+III data set. The 68% and 90% C.L. in- 
terval from the standard two-flavor analysis with equal neu- 
trino and antineutrino oscillations shown by thin and thick 
solid lines, respectively. The 90% C.L. contour for neutrinos 
(least expansive) and antineutrinos (most expansive) from 
the current analysis are dashed. 



trino mixing. To illustrate the difference between neu- 
trino and antineutrino oscillations permitted by the data, 
Fig. [4] shows the allowed regions as a function of the dif- 
ference of the antineutrino and neutrino mixing angles, 
sin 2 2(9— sin 2$, and mass squared splittings, Am 2 — Am 2 . 
Contours have been drawn as in Fig. [2] and a black tri- 
angle near the origin represents the position of the best 
fit. 

In conclusion, a search for evidence of differing neu- 
trino and antineutrino oscillation parameters in the SK- 
I+II+III atmospheric data sample has been carried out. 
The atmospheric mixing parameters for antineutrino os- 
cillations are consistent with those for neutrino oscilla- 
tion. The results agree with the standard SK atmo- 
spheric oscillation analysis, which have also been pre- 
sented. The SK antineutrino oscillation best fit is con- 
sistent with the parameters found by MINOS using a 
predominantly muon neutrino beam. 
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FIG. 4: Allowed differences between neutrino and antineu- 
trino oscillation parameters for the SK-I+II+III data set. 
The 68%, 90%, and 99% allowed region appear in thin, 
medium, and thick lines, respectively. A black triangle 
shows the location of the best fit point from this analysis. 
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